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Abstract 



We present a calculation of the mass of the lowest-lying negative-parity J = 
l/2~ state in quenched QCD. Results are obtained using a non-perturbatively 
C'(a)-improved clover fermion action, and a splitting is found between the 
mass of the nucleon, and its parity partner. The calculation is performed on 
two lattice volumes and at three lattice spacings, enabling a study of both 
finite- volume and finite lattice-spacing uncertainties. A comparison is made 
with results obtained using the unimproved Wilson fermion action. 



I. INTRODUCTION 

The study of the excited nucleon spectrum can provide important clues to the dynamics 
of QCD and the nature of the interactions between its fundamental partons. The observed 
N* spectrum raises many important questions, such as the nature of the Roper resonance, 
and whether the A (1405) is a true three-quark state or a molecular state. For these reasons, 
the study of the spectrum is an important element of the Jefferson Laboratory experimental 
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programme. The phenomenological interest in the excited nucleon spectrum has been com- 
plemented by a flurry of activity in the lattice community. In particular, two calculations of 
the mass of the parity partner of the nucleon have appeared; the flrst employed the highly- 
improved -D234 fermion action whilst the second employed domain- wall fermions [0,0. 
Both calculations exhibited a clear splitting between the masses of the ArV2+ jyl/2- 
states, and the importance of chiral symmetry breaking in obtaining a non-zero mass split- 
ting has been stressed [Q. 

In this paper, we present a calculation of the lowest lying negative parity nucleon us- 
ing an (9(a)-improved Sheikholeslami-Wohlert (SW), or clover, fermion action; preliminary 
results were presented in ref. 0. By choosing the coefficient of the improvement term ap- 
propriately, all 0{a) discretisation uncertainties are removed, ensuring that the continuum 
limit is approached with a rate proportional to a^. 

The calculation of the excited nucleon spectrum places particularly heavy demands on 
lattice spectroscopy. The excited nucleon states are expected to be large; the size of a state 
is expected to double with each increase in orbital angular momentum. Thus a lattice study 
of the excited nucleon spectrum requires large lattice volumes, with correspondingly large 
computational requirements. Furthermore, the states are relatively massive, requiring a fine 
lattice spacing, at least in the temporal direction. These requirements could be satisfled 
with much greater economy using the clover fermion action than using the domain-wall or 
overlap formulation. Thus it is important to establish that the negative parity states are 
indeed accessible to calculations using the clover action. Finally, by comparing the masses 
obtained using the clover action with a calculation, at a single quark mass, using the Wilson 
fermion action, we also gain insight into the nature of the interaction responsible for the 
splitting in the parity doublet. 

The rest of the paper is laid out as follows. In the next section, we introduce the hadronic 
operators and correlators measured in the calculation, describe the fermion action, and pro- 
vide the simulation parameters. Section |T| contains our results for the masses of the lowest 
lying positive- and negative-parity nucleon states using the SW fermion action. Detailed 
discussion, including a comparison with the Wilson fermion case, and our conclusions, are 
presented in Section 0. 



II. CALCULATIONAL DETAILS 

A. Baryon Operators 

For a particle at rest, there are three local interpolating operators that that have an 
overlap with the nucleon: 

A^y^^ = eijk{uJC'y5dj)uk, (1) 

A^2^^^ = eijk{uJCdj)'y5Uk, (2) 

N^^"^^ = eijk{uJC-f4-f5dj)uk. (3) 

The "diquark" part of both A^^i and couples upper spinor components, while that in 
N2 involves the lower components and thus vanishes in the non-relativistic limit M. For 
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some of the lattices in our calculation, the positive-parity nucleon mass is obtained using 
the non-relativistic quark operators [Q, defined by 

^ ^ = i(l + 74)V^, ^"''^ = ^^(1 + 74). (4) 

In practice, lattice calculations confirm the naive expectation that the operators A^i and N3 
have a much greater overlap with the nucleon ground state than N2, and therefore we do 
not use this operator in the fits. 

The correlators constructed from these operators receive contributions from both the 
positive-parity nucleon, and from its (heavier) negative parity partner. However, we can 
achieve some delineation of the states into forward-propagating positive-parity states and 
backward-propagating negative-parity states, or the converse, through the use of the parity 
projection operator (1 ± 74). On a lattice periodic or anti-periodic in time, the resulting 
correlators may be written: 

C^+/-{t) = E ((1 ± 74)a/3(iV.,a(5, t)N,,pm + 

X 

(1 T74)a/3(iV^,a(x,iV, -t)iV,,^(0))) , (5) 

where Nt is the temporal extent of the lattice. At large distances, when t^l and Nt—t ^ 1, 
the correlators behave as 

C^+{t) Ate-<' + Are-'"»^(^*-*) (6) 

i 

Cjv- (t) ^ A-e-'"^* + 4+e-"^^(^*-*) (7) 

i 

where ml and m~ are the hghtest positive- and negative-parity masses respectively in 
channel i. 



B. Fermion Action 

To leading order in a the Symanzik improvement programme amounts to adding the 
well-known Sheikholeslami-Wohlert term to the fermionic Wilson action 

5S = -Cswy E i^{x)(Tf,uFf,ty{x)^jj{x). (8) 

Provided that Cgw is chosen appropriately, spectral quantities such as hadron masses approach 
the continuum limit with a rate proportional to a^. Non-perturbative determinations of Cgw 
have been made in the quenched approximation to QCD in refs. and 

The Sheikholeslami-Wohlert term, eqn. (^, is of magnetic moment form, and it is well 
known that the use of the SW action results in hyperfine splittings that are closer to their 
experimental values than those obtained using the standard Wilson fermion action, see, for 
example, ref. |10|- The SW term also removes the leading chiral-symmetry-breaking effects 



at finite a. In view of these considerations, we compare the splitting between the positive- 
and negative-parity states obtained using the Wilson fermion action at a single light-quark 
mass with the results obtained using the SW fermion action. 
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C. Simulation Details 



The calculation is performed in the quenched approximation to QCD using lattices gen- 
erated by the UKQCD and QCDSF collaborations; calculations of the light hadron spectrum 
using these lattices have appeared in ref. and |[T2| . |r3| respectively. Propagators on the 
UKQCD lattices were computed from both local and fuzzed sources to both local and fuzzed 
sinks; the fuzzing procedure is described in ref. |14]. The parameters used in the calcula- 



tion are listed in Table |. Propagators on the QCDSF lattices were computed using Jacobi 
smearing at both source and sink, described in ref. |]15| . 

The errors on the fitted masses are computed using a bootstrap procedure. In the case of 
the UKQCD data, the same number of configurations are used for each of the quark masses 
for a given f3 = and lattice volume. However, in the case of the QCDSF configurations, 
different numbers of configurations are used at different quark masses even at the same (5 
and volume. This precludes the use of correlated fits in the chiral extrapolations, and thus a 
simple uncorrelated fit is performed, with the uncertainties computed from the variation 
in the x^- 



III. RESULTS 

The masses of the lowest-lying A^^/^+ and N^^'^^ states are obtained from a simultaneous, 
four-parameter fit to the positive- and negative-parity correlators of eqn. (^), constructed 
using fuzzed sources and local sinks (UKQCD) or using smeared sources and smeared sinks 
(QCDSF), using the fit functions of eqns. (|) and (^. We see a clear signal for the mass of the 
negative-parity states, and the quality of a fit is illustrated in Figure |1|; the contamination of 
the negative-parity correlator from the lighter, backward-propagating, positive-parity state 
is clear both in the fits and in the data. The masses of the lightest particle of positive and 
negative parity as a function of on each ensemble are shown in Figures 

For the chiral extrapolations of the hadron masses, we adopt the ansatz 

{amxY = {aMxf + h2{am^f (9) 

where we use upper-case letters to denote masses obtained in the chiral limit, and X is 
either N^/"^^ or N^^'^~ . We include data at different volumes, but at the same /5, in the chiral 
extrapolations, but treat the fuzzed and Jacobi-smeared data independently; we include only 
those masses which are less than of order one in lattice units. The parameters of the fit for 
the positive- and negative-parity states are given in Table ||, and the chiral extrapolation, 
together with the extrapolated masses, shown in Figures 0-^. Note that the fit to m^, rather 
than rux, gives a sensible behaviour in the heavy-quark limit, whilst being formally the same 
at light pseudoscalar masses; a plot of mx vs. at /3 = 6.4 exhibits clear curvature. 

From quenched chiral perturbation theory one expects leading non-analytic terms which 
are linear in m-Tr [T^. The coefficient of the leading term is predicted to be negative. We 



therefore also attempted to fit our data at /5 = 6.4 to the form 

{amxY = {aMxY + bi^am-,,) + 62(a."^7r)^- (10) 
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Fitting the A^^/^+ state we found a positive value of bi, though with a very large error that 
would still accommodate a negative value. Thus it is unclear whether eqn. |TU| provides a 
reliable form with which to extrapolate data obtained with quark masses around that of the 
strange quark. Furthermore, the chirally extrapolated values thus obtained differ from those 
obtained using eqn. ^ by only around 5%. We therefore quote chirally extrapolated masses 
from the fit to eqn. ^ in the remainder of this letter. There has been considerable study of 
the contribution of non-analytic terms arising from pion-induced baryon self energies [jl^ . 
We do not investigate the effect of these contributions in this paper. 

The ideal procedure for obtaining the chirally extrapolated masses would be to apply 
the forms of eqn. ^ and |10| to the infinite-volume limit of the baryon masses obtained at 
each pseudoscalar mass. Unfortunately, our data do not allow this procedure, but we try 
to gain some insight into the possible finite-volume uncertainties by comparing the masses 
otained at a value of the pseudoscalar mass for which we have data on both smaller and 
larger volumes. Specifically, we compare the masses obtained from the Jacobi-smeared data 
at /3 = 6.0 and /3 = 6.2 with a pseudoscalar mass given by 

am„ = 4.8/(L/a), (11) 

where L is the spatial extent of the smaller of the lattices at each /3. The results of this 



analysis are provided in Table |lTl| ; the mass of the negative-parity state on the larger lattices 
is higher than on the smaller lattices at both lattice spacings, by an amount of order 5%. 

In order to look at the discretisation uncertainties in our data, we show in Figure |^ the 
masses in units of rg against the a^/r^, where tq = 0.5 fm is the hadronic scale ||18|,|19|. For 
the lightest positive-parity state, we obtain entirely consistent results at the different lattice 
spacings. For the case of the negative-parity states, there is some trend towards decreasing 
mass at finer lattice spacings. In order to quantify the discretisation uncertainties, we 
perform uncorrelated fits in o^/tq: 



{Mxro){a) = Mxro{a = 0) + ci-) , (12) 

\ro/ 

where X is either iV^/^+ or N^^'^" , yielding 

M^i/2+ro = 2.74(4) 

M^i/2-ro = 4.1(1). (13) 

The fits and extrapolated masses are shown on the figure, together with the experimental 
values of the masses in units of tq. We estimate the difference between the values obtained 
in the continuum limit and that obtained on our finest lattice, f3 = 6.4, as a measure of the 
discretisation uncertainties in our calculation; this is negligible for the positive-parity state, 
and 10% for the negative parity state. 

Though there is a noticeable discrepancy between the lattice and physical values, a direct 
experimental measurement of tq is unavailable, and ro is a better scale for comparing data 
at different lattice spacings than for comparing data with experiment. Therefore we choose 
as our final result the mass ratio of the negative- and positive-parity masses in the quenched 
approximation 
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M^i/2-/Mjvi/2+ = 1.50(3), 



(14) 



where the quoted error is purely statistical, and we estimate systematic uncertainties of 
order 10% due to finite volume effects, and due to chiral extrapolation and discretisation 
uncertainties. This ratio is to be compared with the physical ratio of 1.63. 



IV. DISCUSSION AND CONCLUSIONS 

The calculation exhibits a clear mass splitting between the positive- and negative-parity 
states, in agreement with calculations using the highly- improved and domain-wall fermion 
actions. In view of the dearth of studies of negative-parity baryon masses using the Wilson 
fermion action, it is instructive to compare the mass splitting obtained with the clover 
fermion action with that obtained from the standard Wilson fermion action. A discrepancy 
in the splitting between the two actions would be indicative of lattice artifacts, that would 
presumably be more severe in the case of the Wilson fermion action. In order to study this, 
Wilson quark propagators were computed on the UKQCD 24^ x 48 lattices a.t (3 = 6.2 at 
a quark mass corresponding to irtT^/mp = 0.7; the smearing and fitting procedures were the 
same as those employed in the SW calculation. 

The splitting between the masses of the positive- and negative-parity states obtained 
with the Wilson fermion action, together with that obtained with the SW fermion action 
is shown in Figure The splittings obtained with the two actions are entirely consistent, 
which we will now argue is reasonable. Under the SU (6) spin-fiavour symmetry, the low- 
lying negative-parity baryons, up to around 2GeV, can be assigned to a Z = 1 70-plet. 
Similarly, the low- lying positive-parity baryons can be assigned to an / = 56-plet. Thus 
the splitting in the parity doublets is analogous to the P — S splitting in the meson sector, 
which we know is relatively faithfully reproduced using the Wilson fermion action. Indeed 
an earlier lattice calculation demonstrating that the P-wave baryons are accessible to lattice 



calculation is contained in ref. 20 



Within the / = 1 70-plet, calculations of the masses of these states both in the quark 
model ||2l[], and in large Nq [^, suggest that the spin-orbit contribution is surprisingly 



small, whilst the hyperfine contribution is of normal size, and there is an effective interaction 



carrying the quantum numbers of pion exchange |2^. It is within a multiplet that we might 
expect the choice of action to play a role. 

The continuation of the calculation to lighter quark masses will require careful consider- 
ation. A baryon can change parity through the emission of a vr or r]'; the latter process is 
accessible even in the case of quenched QCD, where the tt and r]' are degenerate in mass, 
as illustrated in Figure The non-unitary behaviour associated with such processes has 



been observed in the scalar correlator p^-p6|, and the presence of such quenched oddities 



in the baryon sector discussed in ref. |T^. Whilst no evidence for non- unitary behaviour 
is observed in this calculation, the lightest quark mass is indeed close to the A^* — > A^vr 
threshold. 



We are grateful to Robert Edwards and Chris Michael for this observation. 
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In this letter it has been shown that both the SW-clover and Wilson fermion actions are 
capable of resolving the splitting between the positive- and negative-parity baryon masses 
in the quenched approximation to QCD, and therefore that these are accessible to relatively 
economical calculation. We obtain a ratio for the masses of the negative- and positive-parity 
states of 1.50(3), where the error is purely statistical, compared with the experimental value 
of 1.63. A more extensive study of the spectrum including the case of non-degenerate quark 
masses, and for the lowest-lying / = § states will appear in a longer paper [E7 . 
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TABLES 



/5 


-T 


L[fm] 




Smearing 


6.4 1.57 


32^ • 48 


1.6 


0.1313, 0.1323, 0,1330, 0.1338, 0.1346, 0.1350 


J 


6.2 1.61 


24^ • 48 


1.6 


0.1346, 0.1351, 0.1353 


F 




24^ • 48 


1.6 


0.1333, 0.1339, 0.1344, 0.1349, 0.1352 


J 




32^ • 64 


2.1 


0.1352, 0.1353, 0.13555 


J 


6.0 1.76 


16^ • 48 


1.5 


0.13344, 0.13417, 0.13455 


F 




16^ • 32 


1.5 


0.1324, 0.1333, 0.1338, 0.1342 


J 




243 • 32 


2.2 


0.1342, 0.1346, 0.1348 


J 


TABLE L The paramet 


ers of the lattices used in the calculation. The labels J and F refer to 


use of Jacobi and "fuzzed" 


quark 


sources respectively. Lattice sizes in physical units 


are quoted 


using ro to set the scale. |n 










6.4 




6.2 6.0 






Jacobi 




Jacobi Fuzzed Jacobi 


Fuzzed 


aM^i/2+ 


0.285(4) 




0.369(3) 0.37(2) 0.512(6) 


0.53(2) 


b2 


2.38(2) 




2.40(2) 2.5(3) 2.42(3) 


2.3(2) 


aMj^i/2- 


0.46(1) 




0.58(1) 0.62(2) 0.89(2) 


0.91(5) 


b2 


2.58(6) 




2.67(7) 2.4(4) 2.1(2) 


1.7(7) 


TABLE IL 


The parameters of the fit of the hghtest positive- and negative-parity states to 


eqn. @. 










P 






amjyi/2- (small) am ^1/2- (large) 


6.0 


0.3 




0.97(1) 


1.01(1) 


6.2 


0.2 




0.664(9) 


0.69(4) 



TABLE in. A comparison of the masses of the lowest-lying negative-parity states on the 
small- volume and large- volume lattices at (3 = 6.0 and f3 = 6.2, at a pseudoscalar mass 
amTT = 4.8/(L/a). 
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FIG. 1. The effective masses for the A^^/^+ channel (circles) and the N'^^'^ channel (diamonds) 
at (3 = 6.2 with k = 0.1351. The lines are from a simultaneous fit to both parities. 



1.5 r 



0.1 




Fuzzed 
Jacobi 

Jacobi: 24^.32 
Jacobi: 16^.32 
O Fuzzed: 16^.48 



0.2 



0.3 



0.4 



FIG. 2. The masses in lattice units of the lowest-lying positive- and negative- parity nucleons 
on the lattices at f3 = 6.0. The curves are from independent fits to the Jacobi-smeared and fuzzed 
data for using eqn. (^), as described in the text. 
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FIG. 3. The masses in lattice units of the lowest-lying positive- and negative-parity nucleons 
at /3 = 6.2. The curves are from independent fits to the Jacobi-smeared and fuzzed data for rrcj^ 
using eqn. (^). 




FIG. 4. The masses in lattice units of the lowest-lying positive- and negative-parity nucleons 
at /? = 6.4. The curves are from fits to using eqn. (|9[). 
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4 ^ 
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1 - 



-a- 



* 



N(939) 



□ Jacobi 
O Fuzzed 



0.01 0.02 0.03 0.04 

(a/r/ 

FIG. 5. The masses of the lowest-lying positive- and negative-parity baryons in units of 
r^^ [010 against a? in units of r^. The lines are linear fits in a^/r^ to the positive- and nega- 
tive-parity b 



O Clover 
O Wilson 



g 0.5 



HI 



m 



0.05 0.1 

(amf 

FIG. 6. The masses in lattice units of the lowest-lying positive- and negative-parity nucleons 
using the SW-clover action at P = 6.2 on the 24^ x 48 lattices (circles). Also shown is the 
corresponding results obtained using the Wilson fermion action (diamonds) on the same ensemble 
of configurations. 
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FIG. 7. Diagram contributing to the decay N'^^'^ — > N^/^^rj' in quenched QCD 



